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ABSTRACT 


Total  strain  energy  of  bcc  to  hep  martensitic  transformation  in 
several  titanium  alloys  has  been  evaluated.  Its  minimum  is  assumed  to 
dictate  the  characteristics  of  a transformation  product  and  various 
features  of  Ti  alloy  martensites  have  been  deduced.  This  analysis  con- 
sidered both  elastic  strain  energy  originated  from  the  transformation 

strain  and  plastic  relaxation  strain  in  relief  of  high  internal  stresses.  ^ 

_ „ \ 

The  primary  mode  of  deformation  evaluated  for  the  latter  was  (101 1 }<101 2> 
twin.  In  Ti-Ta  and  Ti-Zr  alloys,  the  usual  { 334 }g  habit  was  predicted  for 
all  the  twin  variants  and  the  amount  of  twin  shear  was  approximately  0.02. 
The  results  are  similar  to  those  reported  for  Ti  martensite.  In  Ti-Cr  and 
Ti-Mn  alloys,  both  the  { 334}g and  (344}ghabits  were  found  to  be  favorable, 
depending  on  the  specific  variant  of  twin  shear,  the  amount  of  which  was 
approximately  0.04.  Theoretical  predictions  on  the  habit  plane,  shear 
mode  and  twin  fraction  are  in  good  agreement  with  the  experimental  obser- 
vations in  the  literature. 
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1.  INTRODUCTION 


The  importance  of  strain  energy  contribution  to  the  free  energy  change 

during  a phase  transformation  has  been  well  established.  However,  the  explicit 

evaluation  is  often  cumbersome  and  only  a limited  number  of  calculations  has 

been  made.  Over  the  years,  such  simplified  cases  as  a spheroidal  inclusion  with 

hydrostatic  transformation  strain^  have  been  employed  as  approximations. 

Recently,  the  Eshelby  theory  of  the  strain  energy  in  a constrained  transforma- 

tion  has  been  invoked  in  the  treatment  of  more  generalized  transformation  strains 

3 4 5 

involving  martensitic  transformations  of  Ti  , and  steels  ’ . Other  developments 
in  this  field  include  more  elaborate  theories  for  the  evaluation  of  the  strain 
energy  of  an  ellipsoidal  inclusion  having  non-uni formly  distributed  transforma- 
tion strain  in  an  anisotropic  matrix^’'7,  and  the  calculation  of  the  strain 

8-1 0 

energy  in  the  reciprocal  space  . The  latter  has  been  applied  to  studies  of 

11  12 

some  phase  transformation  ’ . 

3 

In  our  first  paper,  we  evaluated  the  strain  energy  due  to  the  bcc  to  hep 
phase  transformation  in  titanium.  The  transformation  strain,  e^,  was  deduced 

1 *5  2 

from  the  Burgers'  shear  mechanism  , and  the  Eshelby  theory  on  the  strain 

energy  of  ellipsoidal  inclusion  was  employed.  It  was  found  that  the  strain 

energy  is  dependent  on  the  orientation  as  well  as  on  the  apect  ratio  of  a 

spheroidal  inclusion.  The  energy  minima  were  found,  when  a disc-shaped  inclusion 

was  lying  on  { 334 } , which  coincided  wvch  the  habit  plane  of  Ti  martensite. 

Subsequently,  the  calculation  was  extended  to  obtain  the  conditions  that 

14 

minimize  total  strain  energy  of  transformation  . It  was  shown  that  even 
when  the  elastic  strain  enerqy  is  m.r.  nzed,  a large  compressive  stress  com- 
ponent within  the  newly  formed  martensite  phase  exists.  This  internal  stress 
can  be  reduced  by  uniform  plastic  deformation  within  the  new  phase.  Denoting 


plastic  deformation  or  plastic  relaxation  strain  (PRS)  as  e^. , total  stress- 

★ 

free  or  eigen  strain,  e..,  becomes, 

* J 


* t r 

e.  . = e. . + e. . 
ij  ij  ij 


(1) 


★ 

Now  the  elastic  strain  energy,  E$,  is  calculated  by  using  e^.,  and  the  work 
done  by  the  plastic  deformation  is  added  to  evaluate  the  total  strain  energy, 
E,  as 


(2) 


where  the  work  done  by  the  plastic  relaxation  strain,  Ep,  can  be  written  as 


Here  is  the  volume  of  a spheroidal  inclusion,  and  tc  and  y are  the  critical 
shear  stress  and  the  magnitude  of  simple  shear  of  the  PRS,  respectively.  The 
basic  hypothesis  of  this  modified  approach  can  be  stated  as  follows:  The 

condition  of  minimum  strain  energy  dictates  the  morphology  and  orientation  of 

a transformation  product  and  the  amount  and  type  of  plastic  relaxation  strain. 

The  theory  was  applied  to  the  study  of  martensitic  transformation  in  pure 
14 

Ti  In  order  to  reduce  the  elastic  strain  energy  resulting  from  the  phase 
change,  a total  of  84  glide  and  twin  systems  in  the  new  phase  was  considered 
as  PRS.  It  was  found  that  for  oblate  spheroidal  martensites  with  aspect  ratios 
(thickness/diameter)  less  than  0.01,  three  types  of  PRS's;  namely,  { T 01 1 }<1 01 2> 
twin  and  <a  + c>  glides  on  (1011 }and  {1122},  produce  strain  energy  minima  at 
shear  strain  of  about  2%.  Theoretical  predictions  on  habit  plane,  shape  shear, 
the  amount  and  type  of  PRS,  etc.  were  in  excellent  agreement  with  experimental 


observations . 


In  this  paper,  we  wish  to  apply  this  modified  approach  to  the  bcc  to  hep 

martensitic  transformation  in  Ti  binary  alloys.  The  experimental  results  of 

habit  planes  for  hep  martensites  have  been  reported  to  be  {334}g,  { 344 } g , or 

both  { 334 } g and  (344}^.  The  most  common  was  the  f 334 > g habit  observed  in  many 

dilute  alloys15,16.  The  { 344} g habit  was  reported  in  Ti-Mn  alloy17,  but  the  mixed 

18  19 

habit  was  also  found  in  Ti-Mn  and  Ti-Mo  alloys.  Because  of  the  availability 
of  experimental  observations,  martensitic  transformations  of  Ti-5.5Cr,  -4.4Mn, 
-7.0Ta,  and  -12Zr  alloys  (alloyinq  content  given  in  atomic  %)  are  selected 
here  for  the  comparison  with  the  predictions  of  the  theory. 

?.  ELASTIC  STRAIN  ENERGY 

Employing  lattice  parameter  data  of  Ti  binary  alloys  collected  by  Hansen 
20  t 

et  al.  , for  all  the  alloys  were  evaluated  again  assuming  the  Burgers' 
shear  mechanism.  The  results  are  tabulated  in  Table  I,  where  the  reference 
coordinate  axes  are  parallel  tc  [TlT]g,  [ZTl]g  and  [01 1 ]g , respectively.  Note 
that  [011]B  is  transformed  to  [0001]H  by  the  Burgers'  shear  mechanism.  Lattice 
parameters  of  a'  phase  in  Ti-Mn  and  -Cr  systems  were  assumed  to  be  identical 
to  those  of  pure  Ti. 

The  main  difference  in  components  of  e!\  among  those  alloys  exists  in 
e^,  which  is  the  normal  strain  component  along  [0001  ]H . In  Ti-Mn  and  -Cr 
alloys,  the  magnitude  of  e^  is  nearly  twice  that  of  pure  Ti.  In  Ti-Zr  and  -Ta 

alloys,  e^.  are  very  similar  to  that  of  pure  Ti. 

Consider  next  the  plastic  relaxation  strain  e^ . In  our  previous  analysis 

of  pure  Ti  martensite,  84  possible  shear  systems  were  considered,  it 
was  found  that  the  shear  systems  can  be  classified  into  four  different  types  in 
terms  of  their  effect  on  E . Two  of  them,  Type  I and  II  were  found  to  lower  E 

a J 

for  thin  plate  morphologies.  For  the  analyses  of  Ti  alloy  martensites,  {1011} 


twin  and  basal  glide  were  chosen  to  represent  Type  I and  II,  respectively. 

The  first  choice  stems  from  the  fact  that  it  has  been  confirmed  experimentally 
to  be  the  primary  martensite  substructure  in  those  alloys.  The  second  is 
included  for  the  sake  of  a comparison. 

Now  for  a given  variant  of  PRS  in  a given  alloy,  the  total  eigen  strain 

★ 

e . . was  obtained  in  the  reference  coordinate  system  as  a function  of  y.  In 

' J 

order  to  determine  the  orientation  dependence  of  E , this  eigen  strain  was 

» I I 

transformed  to  that  in  another  coordinate  system,  Xj,  X^  and  X^.  These  axes 

are  parallel  to  the  principal  axes  of  an  oblate  spheroidal  martensite.  The 

★ ' 

transformed  eigen  strain,  e..  is  given  by 

^ J 


* • * 
e. . = a . a . e 

u ip  jq  pq 


(4) 


where  a.  is  the  direction  cosine  between  the  reference  coordinate  axis  X 
IP  P 

I 

and  the  inclusion  coordinate  Xj.  The  elastic  strain  energy  of  the  oblate 
spheroidal  inclusion  can  be  found  by  the  procedures  described  in  reference  3. 

From  a series  of  computer-aided  calculations,  it  was  possible  to  find  the  disc 
orientation,  which  gives  th  eminimum  elastic  strain  energy,  Em-n,  as  a function 
of  y and  the  aspect  ratio,  k,  of  the  inclusion. 

3.  RESULTS  OF  CALCULATIONS 
3 . 1 Elastic  Strain  Energy  Minimum 

★ 

When  PRS  of  either  Type  I or  II  was  included  in  e..,  the  minimum  elastic 

' J 

strain  energy,  Em^n,  was  reduced  drastically.  Figures  1 and  2 are  the  plots 

E„.  vs.  y for  k = 0 and  k = 0.01,  resoecti vely . The  dotted  lines  are  the  results 
mm 

for  the  Ti-Zr  alloy.  The  Ti-Ta  alloy  behaves  almost  identically  to  the  Ti-Zr 
alloy  for  an  identical  PRS,  reflecting  the  similarity  in  e^..  In  these  two 

* J 

alloys,  the  amount  of  shear  required  for  minimizing  E^  (defined  as  yc)  was 


4 


approximately  0.02  for  Type  I and  0.08  for  type  II,  respectively.  The  solid 

lines  in  Mgs.  1 and  2 show  the  corresponding  variation  in  E . in  the  Ti-Mn 

alloy.  Results  for  the  Ti-Cr  alloy  were  similar  to  those  of  the  Ti-Mn  alloy. 

In  the  Ti-Mn  and  -Cr  alloys,  the  amount  of  shear  needed  to  minimize  . n 

increased  to  0.04  for  Type  I and  0.10  for  Type  II.  When  k = 0,  the  values  of 

were  reduced  to  zero  with  the  addition  of  PRS.  The  lowest  value  of  E . 
mm  mm 

increased  with  a higher  aspect  ratio.  Numerical  results  are  summarized  in 

Table  II.  Note  that  even  at  k > 0,  the  minimum  value  of  E . in  the  Ti-Mn 

mi  n 

alloy  with  Type  I PRS  is  equal  to  that  in  pure  Ti,  which  was  reported  in  re- 
ference 14.  On  the  other  hand,  the  minimum  value  of  E . with  Type  II  PRS  in 

mm 

Ti-Mn  alloy  is  larger  than  that  in  pure  Ti. 

3.2  Total  Strain  Energy  Minimum 

While  E . can  be  reduced  by  a substantial  amount  with  the  addition  of 
mm 

PRS,  this  plastic  relaxation  strain  also  contributes  to  the  total  strain  energy. 
Figure  3 and  4 are  the  plots  of  E vs.  y in  the  Ti-Mn  alloy  with  {1011}  twin 
(Type  I PRS)  at  k = 0 and  k = 0.01,  respectively.  The  critical  shear  stress, 
t , in  Eq.  (3)  for  a given  shear  system  is  expressed  as  cep  where  y is  the 
shear  modulus,  and  the  values  of  a employed  are  shown  in  the  figures.  It  is 
seen  that,  with  Type  I PRS,  the  total  strain  energy,  E,  reaches  a minimum  value 

(Em)  at  y = yc  even  when  a is  as  high  as  0.01  (see  Fig.  3).  Similar  results 
of  E vs.  y were  obtained  with  Type  II  PRS.  In  this  case,  however,  the  total 
strain  energy  was  reduced  only  when  a was  less  than  0.005.  Numerical  results 
of  Em  and  yc  are  given  in  Table  III.  In  the  Ti-Cr  alloy,  the  relationships 
between  E and  y were  similar  to  those  in  the  Ti-Mn  alloy  as  described  above. 

In  the  cases  of  the  Ti-Zr  and  -Ta  alloys,  the  relationships  were  close  to  those 

14 

in  pure  Ti  reported  previously  . 
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4. 


DISCUSSION 


The  total  strain  energy  being  considered  in  our  theory  must  be  smaller 

than  the  magnitude  of  the  chemical  free  energy  change,  |aFc|,  in  the  phase 

transformation.  Here,  lAFjwere  evaluated  by  using  Kaufman's  formul ar  for  Ti 
21 

binary  alloys  . In  all  the  alloys  considered,  it  was  found  that  |AF  |is  greater 

than  E . for  the  inclusion  with  k < 0.01  even  when  PRS  is  iqnored.  Therefore, 
min 

insofar  as  the  energy  balance  is  concerned,  alloy  martensites  with  the  very 
thin  disc  morphology  can  exist  without  any  plastic  relaxation. 

In  order  to  reduce  total  strain  energy  via  the  operation  of  PRS,  it  is 
necessary  that  the  resolved  shear  stress  (RSS)  on  a particular  shear  system 
exceeds  its  critical  shear  stress  (CSS)  at  the  temperature.  The  internal 
stress  due  to  alone  inside  martensites  in  all  the  alloys  considered  was 
found  to  be  essentially  a large  compressive  stress  along  [0001]„.  Let  us  first 
consider  {1011}  twin  system  (Type  I PRS).  The  RSS  on  this  twin  system  under 

the  above  mentioned  condition  can  be  easily  found  by  using  the  transformation 
law  for  second-order  Cartesian  tensors.  For  k = 0,  the  RSS  for  (1011}  twin 

in  a given  alloy  is  given  iri  Table  IV.  These  values  turned  out  to  be  almost 

independent  of  k.  On  the  other  hand,  it  is  difficult  to  determine  the  effect 

of  alloy  addition  on  the  CSS  of  {1011}  twin  system  in  Ti  alloys.  We  estimated 

this  quantity  form  the  plasticity  data  in  polycrystalline  Ti  alloys  at  elevated 

22  23 

temperatures  and  the  CSS  data  in  pure  Ti  * . Table  IV  lists  the  calculated 

24 

values  of  RSS  and  estimated  values  of  CSS  at  temperature  . It  is  seen  that 
in  all  the  alloys,  the  RSS  is  much  larger  than  the  CSS.  Thus,  it  is  expected 

to  observe  this  twin  system  in  Ti  alloy  martensites  in  agreement  with  previous 

. . 17,25,26,27 

studies  ’ . 

The  observed  twin  fraction  in  Ti-4.4Mn  martensite  was  reported  to  be 
0.2  - 0.25^  . Our  prediction  of  the  twin  fraction  is  0.28  ~ 0.29,  when  k is  in 


. 6 - 


the  range  of  0 to  0.01  and  a is  taken  as  0.005.  This  value  of  u corresponds  to 

the  estimated  CSS.  Theory  and  experiment  agree  well  in  this  case.  The  average 

25 

twin  fraction  observed  in  Ti-7.0Ta  martensite  was  0.10  ~ 0.15  . This  agrees 

reasonably  well  with  our  results.  The  predicted  twin  fraction  was  in  the  range 

of  0.17  - 0.20,  corresponding  to  k of  0 to  0.01  and  a of  0.002.  According  to 

the  present  calculation,  martensite  of  the  Ti-Cr  alloy  is  expected  to  exhibit 

the  same  PRS  as  that  of  the  Ti-Mn  alloy.  In  Ti-Cr  alloys,  however,  another 

type  of  twin,  { 1 0T 2 } <T0 1 1 > , was  reported  in  martensite  in  addition  to  the  usual 
27 

{1011}  twin  . It  is  not  clear  why  both  twin  systems  coexist,  since  the  shear 
directions  are  opposite. 

The  RSS  of  basal  glide  (Type  II  PRS)  is  zero  for  k = 0 and  is  always  much 

smaller  than  that  of  (lOTl)  twin  even  for  the  cases  of  finite  aspect  ratios. 

Since  Type  II  PRS  is  less  effective  in  achieving  the  minimization  of  Emin>  the 

operation  of  this  glide  system  is  unlikely  unless  its  CSS  is  considerably 

lower  than  that  of  Type  I PRS.  This  appears  not  to  be  the  case.  Type  II  PRS 

1 7 Pf\  PR 

was  not  reported  in  any  of  the  microstructural  studies 

Finally,  predicted  habit  plane  orientations  are  presented  and  compared 
with  experimental  observations.  The  plane  normal  to  the  minor  axis  of  a thin 
oblate  spheroidal  inclusion  is  defined  as  the  habit  plane  when  the  conditions 
of  the  elastic  strain  energy  minimum  were  reached  by  using  Eq.  (4). 

The  habit  plane  orientations  obtained  without  plastic  relaxation  strain 
were  within  several  degrees  of  [TT3]D.  Those  are  indicated  as  P^  in  Fig.  5. 

In  Ti-Ta  and  -Zr  alloys,  PQ  almost  coincided  with  [4T3]g  and  also  was  essentially 
the  same  as  that  of  pure  Ti.  In  Ti-Mn  and  -Cr  alloys,  PQ  was  shifted  from 
[TJ3]g  toward  [TTl]g  by  3 to  4 degrees.  The  addition  of  plastic  relaxation 
strain  shifts  the  habit  plane  orientation  from  PQ  a way  from  the  001  - 111 
symmetry  line.  The  predicted  orientation  corresponding  to  the  total  strain 
energy  minimum  depended  on  material  (through  e^  and  a)  and  also  on  the  variant 


. 7 - 


of  opera  tint)  IIOTlI  twin.  It  1*.  interest  ing  to  note  that  the  magri  i tude  of  L1" 
was  unchanged  even  though  different  habit  orientations  resulted  with  different 
variants  of  HoTl}  twin  system.  In  Ti-Ta  and  -Zr  alloys,  (10111  twin  as  PRS 
always  resulted  in  the  ( 334 } g habit.  Results  for  Ti-Ta  alloy  are  shown  in 
Fig.  5a  (corresponding  to  Em  at  a = 0,  k = 0).  Six  different  habit  orientations 
resulted  from  the  six  variants  of  this  twin  system  are  indicated  in  the  figure. 

If  a > 0,  Em  was  achieved  at  a smaller  amount  of  y and  the  orientations  predicted 
were  closer  to  PQ  than  the  ones  in  Fig.  5a.  Effects  of  k were  small,  as  the 
predictions  did  not  change  more  than  1 degree  from  the  above  results  with  k = 0 
at  a given  a.  Therefore,  all  the  possible  predicitons  in  Ti-Ta  and  -Zr  alloys 
happened  to  be  the  {334}^  habit  plane,  as  was  the  case  in  pure  Ti14.  No  experi- 
mental determination  of  habit  plane  in  these  alloys  is  available  for  comparison, 
however. 

In  Ti-Mn  and  -Cr  alloys,  two  out  of  the  six  variants  of  {10Tl}  twin  also 
produced  basically  the  [433]^  habit.  In  Fig.  5b,  the  results  for  Ti-Mn  alloy 
are  shown.  Two  of  the  six  twin  variants  produced  the  habit  orientations  within 

B‘ 

The  specific  Kl  ^P^nes  of  these  variants  are  indicated  in  Fig.  5.  These  four 

orientations  can’ be  classified  as  the  { 344}g  habit.  These  were  determined  at 

the  total  strain  energy  minima  with  a - 0 and  k = 0.  As  mentioned  previously, 

k hard^^fected  the  habit  plane  orientation,  but  effect  of  a was  significant. 

When  i was  increased  to  0.005,  shifts  in  the  habit  orientations  toward  P was 

o 

observed  as  s..own  in  Fig.  5c.  Two  variants  of  { 1 OT 1 } twin  still  produced  the 
{ 344 } g habit.  Results  for  Ti-Cr  alloy  were  essentially  identical.  Thus  in 
these  alloys,  two  groups  of  habit  orientations  are  anticipated;  one  is  close  to 
the  f334}„  habit  and  the  other,  the  (344id  habit.  This  is  in  agreement  with  an 
experimental  finding  in  Ti-Mn  alloys.  Liu  and  Margolin  reported  these  two 

1 g 

habits  in  alloys  containing  4.35  to  5.25  wt.Mn  . However,  in  Ti-Cr  alloys, 


a few  degrees  of  The  other  two  habit  orientations  were  close  to  [434] 


- 8 - 


i n disagreement.  with  the 


only  I 334 1 B habit  was  reported  experimentally 
above  prediction.  In  this  alloy,  a second  type  of  twin,  ;10T2}<101 1--,  was  also 
observed  in  the  martensite*^.  The  origin  of  the  second  shear  is  not  known, 
but  our  calculations  suggest  that  lowering  of  the  strain  energy  is  unlikely  to 
be  the  reason. 
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TABLE  II.  Minimum  values  of  Em1- n and  corresponding  amounts  of  PRS 
required  to  achieve  the  strain  energy  minimization,  y . 


Alloy 

i 

R S 

— 

S 

Estimated  CSS  at  M$ 

T 

a 

Tc 

a 

Ti-4.4  Mn 

MPa 

687 

0.0178 

MPa 

157-470 

0.0041-0.0122 

-5.5  Cr 

627 

0.0169 

124-353 

0.0034-0.0095 

-12  Zr 

304 

0.0098 

71 

0.0023 

-7.0  Ta 

284 

0.0086 

74 

0.0023 

Pure  Ti 

294 

0.0099 

50 

0.0017 

Table  IV. 


Resolved  Shear  Stress  (RSS)  due  to  Transformation^  Strain 
and  Estimated  Critical  Shear  Stress  (CSS)  of  (1011 }<T01 2> 
Twin  in  Ti  Alloy  Martensites. 
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